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Abstract 

The 5N formula that relates the final curvature perturbation on comoving slices to the 
inflaton perturbation on flat slices after horizon crossing is a powerful and intuitive tool to 
compute the curvature perturbation spectrum from inflation. However, it is customarily 
assumed further that the conventional slow-roll condition is satisfied, and satisfied by all 
components, during horizon crossing. In this paper, we develop a new 5N formalism for 
multi-component inflation that can be applied in the most general situations. This allows us 
to generalize the idea of general slow-roll inflation to the multi-component case, in particular 
only applying the general slow-roll condition to the relevant component. We compute the 
power spectrum of the curvature perturbation in multi-component general slow-roll inflation, 
and find that under quite general conditions it is invertible. 



1 On sabbatical leave from Department of Physics, KAIST, Daejeon, Republic of Korea 



I. INTRODUCTION 



Thanks to recent advances in observational technologies, an era of precision cosmology has 
commenced, and there have been many studies attempting to determine, or at least constrain, 
possible models/theories of the early universe by cosmological observations. Notably, the 
analysis of the first year WMAP data strongly indicates that the universe is spatially flat, 
that the primordial perturbation is adiabatic and Gaussian, and that the spectrum is almost 
scale-invariant [1]. By and large, these support the now-standard picture that our universe 
experienced an inflationary phase at its early stage. Furthermore, the PLANCK satellite, to 
be launched within a couple of years from now, will provide us with more detailed information 
about not only the CMB temperature spectrum but also the polarization spectrum [2]. 

On the theoretical side, both supersymmetric particle theory and string theory lead to the 
expectation that there are many scalar fields. It is then important to understand the cosmo- 
logical dynamics of mult i- component scalar fields and formulate a theory of the cosmological 
perturbations generated from vacuum fluctuations in those fields during inflation. 

The inflationary dynamics of multi-component scalar fields has been investigated by var- 
ious authors from various aspects. In particular, the SN formalism [3-5], in which the final 
curvature perturbation on comoving slices is expressed as the e-folding number perturba- 
tion, which in turn is given in terms of the inflaton perturbation on flat slices after horizon 
crossing, has proved to be a powerful tool to compute the curvature perturbation spectrum 
from inflation. 

However, most of the previous work [3, 6, 7] has been based on the assumption that all 
the fields satisfy the conventional slow-roll condition. The conventional slow-roll condition 
assumes that the slow-roll parameters are not only small but also slowly varying. In the 
case of a single inflaton field, the slow-roll parameters should be small in order to produce 
a scale-invariant spectrum but there can be cases in which the slow-roll parameters are not 
slowly varying. In such cases, the resulting power spectrum of the curvature perturbation 
can be quite different from the conventional slow-roll spectrum. Of course, since the slow- 
roll parameters are small, the spectrum is still almost scale-invariant so that it is consistent 
with present observational data. Thus it is important to relax the conventional slow-roll 
condition and investigate the power spectrum of the curvature perturbation generated from 
general slow-roll inflation [8-11], which in the case of a single inflaton field means assuming 
the slow-roll parameters are small but not necessarily slowly varying. 

In the case of a multi-component inflaton, one has many slow-roll parameters but only the 
relevant ones need to be small to ensure an approximately scale-invariant spectrum, i.e. one 
should only apply the general slow-roll condition to the relevant component, with essentially 
no constraints on the others. This makes the task of calculating a general formula for the 
spectrum challenging. 

In this paper, we formulate the cosmological perturbation of a universe dominated by a 
multi-component scalar field during inflation in a form as general as possible, extend the SN 
formalism developed previously [3] by introducing a new quantity SN, and derive a formula 
for the power spectrum of the curvature perturbation in multi-component general slow-roll 
inflation. 

This paper is organized as follows. In Section II, we briefly review the SN formalism [3], 
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on which our analysis relies heavily. That is, we introduce the perturbation in the e-folding 
number, SAf, and show that 5Af between two different perturbed spacelike hypersurfaces is 
equal to the difference of the curvature perturbations on the two hypersurfaces on super- 
horizon scales. Taking the initial hypersurface to be flat and the final one to be comoving, SAf 
becomes equal to the curvature perturbation on the final comoving hypersurface, TZ C , which is 
conserved in time on scales sufficiently greater than the horizon size if the final hypersurface 
is taken at a sufficiently late time that the perturbations have become adiabatic. 

In Section III, we first write down the perturbation equation for a multi-component scalar 
field in the flat slicing. We introduce a quantity SN which closely resembles the perturbation 
of the background e-folding number, N, which we regard as a function of <fi and <fi. We then 
derive an evolution equation for SN valid from sub-horizon scales to super-horizon scales. 
Finally, we show that this SN can be identified with SAf on super-horizon scales, hence with 
the curvature perturbation on the comoving slice, 7Z C , at late times. 

In Section IV, using the equations derived in the previous section, we derive the general 
formula for the power spectrum of 1Z C in multi-field general slow-roll inflation to leading 
order. In Section V, we conclude our paper. 

Some mathematical formulae used for manipulations of the multi-component perturbation 
equations and some connections with previous formulae are given in the Appendices. 



II. SAf 

Following Ref. [3], we introduce the perturbation of the e-folding number SAf and show 
its relation to the curvature perturbation. 
The background metric is 

ds 2 = dt 2 - a(t) 2 5ij dx*dx j . (2.1) 
We define the e-folding number as 

N = I Hdt , (2.2) 

J *fin 

where H = a/a and t^ n is a late time when all trajectories have converged, i.e. a time after 
complete reheating when the curvature perturbation TZ C has become constant; see Eq. (2.11) 
below. 

We write the scalar part of the perturbed metric as [12] 

ds 2 = (1 + 2A)dt 2 - 2diB dt dx i - a 2 [(1 + 271)5^ + 2a~ 2 d i d j E] dx i dx j . (2.3) 
The perturbed e-folding number is defined by 

N = f \edr, (2.4) 

where r is the proper time, dr = (1 + A)dt, and 9 is the volume expansion rate of the 
constant time hypersurfaces. We write 6 as 

\e = H(l + K). (2.5) 
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In the rest of this paper, we expand the perturbation variables in terms of spatial scalar 
harmonics Y k (x l ) with the eigenvalue k 2 [12], e.g. A(t,x l ) = J2 k A k (t)Y k (x l ), and deal with 
the harmonic coefficients by omitting the eigenvalue indices. For example, we denote A k 
simply by A. 

The perturbation of the expansion rate /C is expressed in terms of the metric perturbations 
in Eq. (2.3) as 

}C = -A + ±<R,-^S , (2.6) 

where q 2 = k 2 /a 2 , 

S = E- 2HE - B , (2.7) 

and 1Z describes the intrinsic curvature perturbation of the constant time hypersurfaces, and 
S the shear of the vector unit normal to them. 

We note that, assuming the matter anisotropic stress is negligible, S satisfies the equation, 

S + HS = A + TZ. (2.8) 

Thus, apart from the modes of S sourced by A + 1Z that remain regular in the limit q 2 — > 0, 
q 2 S decays rapidly as a~ 3 once it is outside the horizon, and this decaying behavior is 
independent of the choice of gauge. 

The perturbation in the e-folding number is written as 

5X(t &n , t) = M - N = K{t) - TZ(t &n ) -\f q 2 S dt . (2.9) 

Taking the initial, t = t ini , hypersurface to be flat and the final, t = tfi n , hypersurface to be 
comoving 

<W(* fin , t ini ) = -ftc(*fin) ~ \ [ ^ <?S dt . (2.10) 

As discussed in the previous paragraph, q 2 S is either rapidly decaying or S is regular in the 
limit q 2 — > 0. Hence, the second term in the right hand side of this equation is negligible 
on scales sufficiently greater than the horizon. So, taking the initial time t ini sufficiently late 
that this term is negligible, we have 

&/V(tfin,fini) ^ -ftc(ffin) ■ (2-H) 

This relation is valid irrespective of whether the background universe is dominated by scalar 
field or not. 



III. MULTI-COMPONENT SCALAR FIELD DURING INFLATION 

In the inflationary universe, the curvature perturbation is generated from vacuum fluc- 
tuations in the scalar field when the wavelengths of the Fourier modes are well inside the 
horizon. Therefore generically it is necessary to solve the evolution of these vacuum fluctu- 
ations until they are well outside the horizon. However, for multi-component inflation, the 
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curvature perturbation 1Z C can still vary in time on super-horizon scales until the background 
trajectories in the phase space converge to a single trajectory. The convergence may occur 
only after the universe becomes totally radiation-dominated, after which the background 
universe undergoes a universal evolution and the only remaining modes of scalar-type per- 
turbations are adiabatic ones (we do not consider the case of isocurvature perturbations that 
persist until the present). It is at this stage that 1Z C becomes constant in time. 

Thinking of the fact that what we need is only the final value of the curvature pertur- 
bation, 7£ c (tfi n ), it seems highly redundant to know the evolution of all the components. In 
other words, if we can identify the part of the perturbations that contributes to 7£ c (£fi n ), we 
may solve only that part to obtain its value without solving the full perturbation equations. 

What the 8J\f formula (2.11) tells us is that this final amplitude of the adiabatic pertur- 
bation 7£ c (ifi n ) is given by 5Af(t ini , £ fin ), i.e. the perturbation in the e-folding number between 
the initial flat hypersurface at t — £ ini during inflation and the final comoving hypersurface 
at t — tfin when TZ C has relaxed to a constant. Thus, the essential information we need is 5J\f 
on super-horizon scales during inflation. Therefore, it will be convenient if we can introduce 
a quantity that can be evolved from sub-horizon scales to super-horizon scales and that can 
be identified with 5Af there. The purpose of this section is to introduce such a quantity 
which we denote by 5N in a multi-component scalar field dominated stage, and argue with 
several pieces of strong, convincing evidence that it indeed matches to 5Af on super-horizon 
scales. A rigorous proof of, including any necessary conditions for, the equivalence between 
8M and SN on super-horizon scales will be a topic of a subsequent paper [13]. 

The reason why we denote this quantity by 5N is that its definition reduces to the 
perturbation of the background e-folding number in the super- horizon limit. 



A. The background 

In this subsection, we review the evolution equations for the background, and work out 
the properties of the e-folding number A as a function on the phase space of the scalar field. 
We denote the dimensions of the field space by D. 



1. Basic background equations 



We assume that for t < t ini , i.e. while modes are leaving the horizon during inflation, we 
can take the action to be 




(3.1) 



We do not require this effective description to continue to be valid for t > t m i. Indeed, we 
expect that in most cases it will break down some time before tg n . 

Here /i a b is the metric of the D-dimensional field space, and following Wald [14], we adopt 
the abstract index notation, i.e. we use boldface Latin indices a, b, ... to denote tensors and 
their contractions in the field space, while the first few Greek letters a, (3, . . . will be used 
to label the components of tensors with respect to a basis of vectors e^ and covectors e^ 
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in the field space, hence a = 1,2, . . . , D. Practically, one may regard the abstract index 
notation as the component notation for an arbitrary basis. One advantage of the abstract 
index notation is that it makes expressions involving derivatives in the curved field space 
concise. Thus, for example, if we take Q be the coordinates on the field space and use the 
usual mathematical definition of vectors as derivatives, we have v a <-> v a d/d(j) a where d/d<p a 
is the coordinate basis. 

The homogeneous background equations are 

a + 3H<j) a + h ah V h V(<f)) = , (3.2) 
SH 2 = i W a b + V , (3.3) 

H = -^ ab a b , (3.4) 

where V a is the covariant derivative on the field space, a = D<p a /dt, and D/dt is the 
covariant time derivative along a background trajectory in the field space. 

Because the background trajectories are determined by specifying both <p a and <p a in 
the field space, it is convenient to extend the covariant derivatives to the phase space in 
which <f) a and <f) a are regarded as independent coordinates. The covariant derivatives in the 
phase space are defined in Appendix A. We denote them by the subscripts <p a and <p a , for 
example Qp and Qp, and they are used extensively below. However, readers who are not 
particularly interested in the curved field space generalization may just regard them as the 
standard partial derivatives with a, b, . . . being interpreted as the component indices. 



2. Np and 

For t < ii n i, we define that iV is represented in phase space (0, 0) as 

rt{4>,4>) 

N((f),<j)) = Hdt , (3.5) 

where the integral is performed along the trajectory that passes through ((f), (p). Therefore, 
we can write the Hubble parameter as 

A> b + N^ h = H , (3.6) 

where the covariant partial derivatives are defined in Appendix A. Differentiating this equa- 
tion with respect to a and a respectively, we obtain 



DN, 



= -N^ (>) a - A^ bJ R b cda c d + Hp 



dt 

= Ay/i b %c a + 3N^ h H^ - A 0b i? b cda c d + Hp , (3.7) 
= -Np + 3HNp + ZN^ h Hp + Hi & , (3.8) 
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where -R b c da is the scalar field space curvature tensor. From Eq. (3.3), we have 
Using these equations, we obtain the evolution equations for N^ b and N^b as 

D 2 ( h ab N \ D / h ab N \ f h de N ; e \ . fh cd N^ c 



dt 2 a 3 dt\ a 3 Dca " M a 3 * M a 
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1 D /a 3 a b \ , fh cd Ni d \ 1 D 



/l ab ^b D f h ah N; b \ iV; b b a > 

+ (3-11) 



a 3 dt \ H ] \ a 3 / 3a 3 dt\H 

and 



a 3 V a 3 / 2#a 3 6#a 3 

Here, Eq. (3.10) has the particular solution 

a 3 ~~ QH J a 3 ^ 12j 
The following contraction, derived from Eq. (3.8) using Eq. (3.6), will also be useful 



B. Perturbations 



In this subsection, we review the perturbation equations for a universe dominated by a 
multi-component scalar field, and introduce a quantity 5N which we will identify with 8H on 
super-horizon scales and whose sub-horizon to super-horizon evolution is particularly simple. 

1. 5 (ft on flat slices 



During inflation, the metric perturbations A and S satisfy 



A = -h ab ^-5(f) D + 
Z ti 



K 

H ' 



q 2 S 



\h ah 



D 

Jt \ H 



H ■ q* 



(3.14) 
(3.15) 



We will use a subscript f to denote a quantity evaluated on flat hypersurfaces, i.e. in the 
gauge 1Z = 0. In particular, 



50? 



5<p a - , 



(3.16) 
(3.17) 



These variables are gauge invariant. 

The scalar field perturbations on flat hypersurfaces satisfy [3] 

S'jf + SHSjf - R\ cd <p h <p c 5tf + q 2 S<Pf + h^V^Stf = ~ (^-f^j h hc 6tf . (3.18) 
In the limit q 2 — * 0, Eq. (3.18) has the solution 

Stfoc^, (3.19) 

corresponding to the super-horizon adiabatic growing mode. Comparing with Eqs. (3.10) 
and (3.12), we see that, in the limit q 2 — > 0, Eq. (3.18) also has the solution 

h ah N, b 6 a f l dt 



5cf>f oc --?-/-, (3.20) 

and, using Eq. (3.11), 

• h ah N^ b IDf^f'dt 

Taking the lower bound of the time integral, t*, in the above equations to be a sufficiently 
late time, this solution describes a super-horizon decaying mode. Note that a different choice 
of t* corresponds to adding the growing mode solution (3.19) to it. Note also that the time 
integral is defined even after inflation when the universe is no longer dominated by the 
scalar field. Thus as a simplest choice we set t* = t &11 , the time at which we evaluate the 
final amplitude of the curvature perturbation 1Z C . 



2. Wronskian 

The Wronskian, W, of solutions 54>i,i and 64>f,2 of Eq. (3.18) is defined as 

W(S^ 2 , 50 fil ) = h ah (<% 5^ - 8<i>l 2 6$^ , (3.22) 

and has the property 

W oc \ . (3.23) 

cr 

We also define the super-horizon Wronskian, Wo(50f,o, as the Wronskian of a solution 
5(f)f of Eq. (3.18) and a super-horizon solution S4>i,o of Eq. (3.18), i.e. a solution of Eq. (3.18) 
in the limit q 2 — > 0. Then 

j t {a 3 W ) = q 2 (a 3 /i ab <^ ) 5$ (3.24) 

and 

{a 3 W ) + 5Hj t {a 3 W ) + q 2 {a 3 W ) = 2q 2 j f (a 3 /i ab <J0f >o ) 5$ . (3.25) 
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These evolution equations for a 3 W motivate our final choice for 5N f , Eq. (3.39), introduced 
in the next subsection. 

Taking the super-horizon decaying mode of Eqs. (3.20) and (3.21), 



h ah N! 



+ 



QH 



dt 



*fin 



and 



8ft 



f,d 



/i ab iV> N^ft°ft a id (ft s 



a 3 2Ha 3 ' 6 dt \ H 

we get the super-horizon Wronskian 

a 3 W (5ft { , d , S(f>{) 
= N^ft? + N^Sftf - ^N^ft a ft h 5ft* + h ah 

= N^5ft? + iV> [d'ftf - 0M f ) + l -q 2 [a 3 J* ^ St . 



1 dt 



*fin 



dt\H Vi H 



(3.26) 



(3.27) 



(3.28) 



3. 5N 

In this section we introduce a quantity 5iV which has the key properties that it behaves 
as 8M on super-horizon scales, but is locally defined apart from the background function N 
and so is straightforward to evaluate. 2 Although not mandatory, we also look for 5N whose 
sub-horizon to super-horizon evolution becomes as simple as possible. 

With this motivation, we define 3 



5N = fV50 a + N^0 a 

2 



where 



and 



X D 



2H * T ra 18// eft \ H 

X (hp. 



18//// ' 

with X being some, as yet unspecified, function of t. 



(3.29) 
(3.30) 

(3.31) 

(3.32) 



2 Strictly speaking, our final choice for X in Eq. (3.38) is not local, but it only depends on background 
quantities and is well approximated by a local quantity. 

3 See Appendix D for the relationship between this definition and the configuration space formula used in 
Ref. [3]. 
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There are several possibilities for the choice of X. The choice X = gives 



5N = N^5(f) a + (<50 a - (3.33) 

and 

IV = iV> . (3.34) 

Note that we may interpret 



5<f) a - (j) a A = 5 j . (3.35) 

Hence 5N may be regarded as the perturbation of the background e-folding number, pro- 
vided that {5(j) a , 5(d(j) a / dr)} is regarded as the perturbation of the background field, i.e. the 
difference between two adjacent background trajectories in the phase space. This choice 
corresponds to the simplest extrapolation from the background SN with (0, <fi) as the phase 
space variables. Note also that, from Eqs. (A20) and (A24), one has a similar interpreta- 
tion if one chooses (0, ir) as the phase space variables, where 7r a = <fi a /H, provided that 
{5(f) a , 5(d<p a / dM)} is regarded as the perturbation of the background field. 

Another choice may be motivated by taking (0, 7r) as the phase space variables with 
5n a = 5(f) a /H. Eq. (A19) then suggests the choice X = 3N wa iT a . This gives 

SN = N^5(f) a + N n * \^n a - 7r a ^ ) (3.36) 

and 

N^a = . (3.37) 

Our final choice is 

X = - 3Ha 'J tti fz =1 + [lp) ' (3 - 38) 

relevant in the case of flat slicing, and motivated by Eq. (3.28) for the super-horizon Wron- 
skian a 3 W . With this choice of X, we have 

5N = N r 5<p a + AV (5'<P a - 0M) + p 2 (a 3 J* ^ St , (3.39) 

= N r~\ (« 3 f d i) ^4 (3 - 40) 




*fin 

and 

N^a D /Nj 



a 3 dt V a 3 



(3.41) 



The choice of X is irrelevant on super-horizon scales, as long as X does not grow rapidly 
in time, but it affects the evolution of 5N from sub-horizon to super-horizon scales. We will 
use this freedom to optimize this evolution, but for the moment let us leave the function X 
in 5N unspecified. 
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Choosing flat hypersurfaces, we have 



SN f = 5N-TZ 



(3.42) 



Using the gauge transformation property of SN, one can readily show that this is a gauge- 
invariant quantity. Taking the time derivative of SNf, we find 



SN f 



3 q 



(3.43) 



where e = —H/H 2 is one of the slow-roll parameters. Thus SNf is constant on super-horizon 
scales. In the special case that our scalar field description is still valid at t = t fin , i.e. if it is 
valid beyond the point where convergence of trajectories occurs, we can evaluate the constant 
by inserting the super- horizon adiabatic growing mode, Eq. (3.19), into the definition of S N, 
Eq. (3.30), to get 



SN { (t ini ) = SNf{t &n ) 



H<p s 



— 7^c(£fin) 



(3.44) 



*=*fin 



Comparing this with SM given in Eq. (2.11), we see that (57Vf (t ini ) = M/"(tfin, £mi), where t m \ 
is a time when the mode is sufficiently far outside the horizon. Thus we have established 
the equivalence of SNf and SM on super-horizon scales in this special case, and the natural 
form of Eqs. (3.33) and (3.36) and their equivalence convinces us that this is true generally. 
We aim to prove this rigorously in a subsequent paper [13]. 

Taking the derivative again we get the equation of motion for SNf 



SNf + 5HSNf + q 2 SNf = -2q 2 N^c 
I I I - — I X 



-Hq 2 
3 y 



l- 6 - 
3 



e 
QH 



1 -2e 
QH 



X + 



QH 2 



X 



(3.45) 



Now let us fix the function X. As we may expect from the evolution equations for the 
super-horizon Wronskian a 3 W , Eqs. (3.24) and (3.25), if we choose X as in Eq. (3.38), the 
Sf terms in Eqs. (3.43) and (3.45) disappear: 



SNf = -q^Mf 



(3.46) 



and 



SNf + 5H SNf + q 2 SNf = -2g 2 N^ a ^ . (3.47) 

What we want to do is to solve the above equation from well inside the horizon to well 
outside the horizon, in particular to the time t = t ini at which SNf and SM can be identified 
with each other. The left hand side of Eq. (3.47) gives a simple evolution equation for the 
single quantity we are interested in, but the right hand side is not so simple. In general, 



we would need to know the solution for 



However, then we would not need to solve 



Eq. (3.47) in the first place, because SNf can be readily constructed from S<pf . Therefore, 
we look for situations in which the whole knowledge of S(f)f may not be needed. For this 
purpose, we decompose S(pf into relevant and irrelevant components, 



(3.48) 
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with the relevant direction defined by N^ a , 

S^^P^Stf, 5<t>l = Ql5<t>\ (3.49) 

where 

so that, from Eq. (3.46), 5(j)* can be expressed entirely in terms of SNf. Then Eq. (3.47) can 
be rewritten as 



SNf + 5H 5N t + q 2 5N { = 2 ( j^fjf; ) 6N t - 2q%J ( j>l . (3.51) 



Thus, in situations when the contribution of the S(j)\ term can be neglected, the equation 
forms a closed system, making it possible to solve for SNf without the full knowledge of the 
scalar field perturbations. Furthermore, we can use observational constraints on the spectrum 
to force SNf to obey the general slow-roll condition, while we have no real constraints on 
S(f>± apart from that its contribution should be small. 

Before closing this section, let us further rewrite Eq. (3.51) in a form suitable for explicit 
calculations. Changing time variable to x — k£, where £ = — J dt/a is minus the conformal 
time, we obtain 

SN'J --5N'+5Nf = -^- SNl - 2N ; a 50* (3.52) 

x x n v 

where 

^"(^e) 2 ^^' (3 - 53) 

SNf = dSNf/dx, II' = dll/dln£, and a£ ~ for e <C 1. This form of the equation for 
SNf will be used in the next section. 



IV. MULTI-FIELD GENERAL SLOW-ROLL POWER SPECTRUM 

At sufficiently early times when all cosmologically relevant wavelengths are well inside the 
horizon scale, we assume that the time scale for the background scalar field is of the order 
of the Hubble time and there exists a set of orthonormal basis vectors, e\ (i = 1, 2, • • • , D), 
in the field space that remain approximately constant in time. In other words, we assume 
that all the frequencies associated with quantum fluctuations of the scalar field are much 
larger than H and their components with respect to the basis e\ can be regarded as mutually 
independent. 

Under the above assumption, we quantize the scalar fields as 

D 

<tyf(M) = [ai(k)<W(M) + <4(-k)(ty? t (M)] , (4-1) 
i=i 
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where Oj(k) and aj(k) are annihilation and creation operators, and the scalar field mode 
functions 8(f)f(k,t) are normalized as 

(4.2) 
(4.3) 



h ah \5<f)t(k, t) 6$(k, t) - 5tf{k, t) 6<%(k, t)\ = , 
K b ^6(f)f\k, t) 6<%(h, t) - 5(f>f(k, t) S^(k, t)] = 2iq 3 5 ij 
The mode functions have the asymptotic behavior 

lim 6$ = qe lx et, 

X—fOO 

lim 8(f>f = -iq 2 e ix ef, 

x^oo 

with 

D 

h ah e*e) = 5 l3) ^^ = h' 



,ab 



Then the creation and annihilation operators obey the quantization conditions 

1 



(4.4) 
(4.5) 

(4.6) 
(4.7) 



[a i (k),a i (l)] = 



aj(k),a](l) 



(4.8) 



Given the quantized scalar field perturbations S(f>f, the quantized 5N{ is expressed as 



D 



5N{(k,x) = M k ) SN i (k,x) + al(-k) 5N*(k,x) 



i=i 



where, from Eq. (3.29), SNi(k,x) is given by 
Eq. (3.52) translates to 



5N » _ t SN > + 6Ni = —SN'i - 2%.50f ± , 
x x 11 v 



with the asymptotic behavior 



lim 6Ni = -iq 2 e ix ni^ 



(4.9) 

(4.10) 

(4.11) 
(4.12) 



following from Eqs. (4.4) and (4.5). 
The power spectrum is given by 

In 2 i 
~ -P(k) 5{k - 1) = lim (6N{{k, x) 5N\{\ x) 



k 3 



x^O 



Therefore 



v ' i=i 



(4.13) 



(4.14) 



Thus what we have to do is to solve Eq. (4.11) for SNi with the asymptotic boundary 
condition (4.12). 
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A. Zeroth order 

To zeroth order, we neglect the terms on the right hand side of Eq. (4.11), as they generate 
deviations from scale invariance, leaving 



SNf - UNf + <wf = . 



x 



This equation has mode functions 



1 



n (x) — ( 1 - ix - ^x 2 ) e lx . 



(4.15) 



(4.16) 



Properties of n (x) are given in Appendix B. The solution for with the asymptotic 
behavior (4.12) is 

5N ( i °\k, x) = 27ri II^oo) n (z) , (4.17) 

where 



Note that 



IT EE 



D 



2vra 2 £ 2 



(4.18) 



i=i 

where II 2 was defined in Eq. (3.53). Therefore the zeroth order power spectrum is 

P{k) = n 2 (oo) . (4.20) 

Note that at this order II 2 is regarded as constant and so the evaluation at x — ^ oo can be 
dropped. 



B. First order 



The Green's function solution of Eq. (4.11) is 
5Ni(k,x) = 5Nf\k, x) + 18 / ^ {Re [n (x)\ Im [n (u)\ - Im [n (x)\ Re [n (u)}} 

Jx u 



X 



irr 



(4.21) 



Taking the limit x — > and using Eq. (BIO), we get 

\mi8NAk,x) = lim5Af } (A;,x) + 181im / — Im[n («)l x 5 Re[n (u)] 

x^o x^o x^o J x u y 45 



X 



(4.22) 
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fl. 



As long as x 2 U' /U and a; 2 N^ a 5(f)f ± become negligible as x — > before IT 7 /]! and N^ a 
stop being small, i.e. before our approximation scheme breaks down, then we can complete 
the limit x — > without needing to apply distinct superhorizon methods [13, 15]. Assuming 
this and completing the limit we get 



7T Im M«)] 

Hi 



irr 

u TT 



(4.23) 



Substituting this into Eq. (4.14), we obtain 

D 



i=i 



5Nf > \k,0) + 18 I ^Im[noH] 



du 



if 



~5N[(k,u)--S^6<f,f L (k,u) 



(4.24) 

Keeping all the potentially leading order corrections from each term on the right hand side 
of Eq. (4.11), we have 

D i 2 
{2ir) 2 P{k) ~ J] \SNf\k, 0) 



i=l 



+ 18j2^\k,0) r^ImKH] 
i=i Jo u 

D 



I^AT< ( "'(fc,u)-N^« L (*,«) 

w J. J. 



+ C.C. 



— lm[n (u)}NiJ^f ± (k,u) 



Defining 



and 



z\\(x) 
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Im [n (a;)] Re [n()(x)] 



2±(a0 = — Im i n o(x)} 

Ob 

with properties given in Appendix B, we have 

p(fc) = n 2 (oo) + 2n 2 (oo)y o ^ -(0 

1 D 



(4.25) 

(4.26) 
(4.27) 



(4.28) 



For scale invariance of the spectrum, we require II 2 to be approximately constant. A 
reasonable extrapolation of this is to assume that N; a is approximately constant, and hence 

N^ a is small. This would naively make the third term in Eq. (4.28) first order and the last 
term second order. However, from Eqs. (3.49) and (4.18), 



D 



D 



En 



i=l 



2ira 2 £ 



(4.29) 



i=i 
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Using Eqs. (4.4) and (4.6), we see that this vanishes in the limit £ — > oo. If <50| and 5<p\ only 
weakly mix, then the third term in Eq. (4.28) will be further suppressed making it naively 
second order. Then the last two terms in Eq. (4.28) could both be regarded as second order, 
though as now at first order we have no direct constraints on <50^ one should be somewhat 
cautious about this. 

If we neglect the 64>j_ terms, Eq. (4.28) reduces to 



\nP(k) = lnn 2 (oo) + 2 J -± z\\{k£) -(£) 
= lnn 2 (oo) + 



Jo 



■KW(kO] 



in it 



2W 
+ 3 11 



(0 



(4.30) 
(4.31) 
(4.32) 



where the definition and properties of W(x) are given in Appendix B. This has a form 
equivalent to the single component case [8-10], though with a different interpretation of the 
inflationary source In II 2 + 2IT/3I1. 4 



C. Inverse formula 

In the special case when the spectrum has the form of Eq. (4.32), it is possible to invert 
the spectrum to obtain the function II 2 , in the same way as was done in the single field 
case [11]. 

Let us introduce the function 



s{x) = - 

7T 



3 3cos(2x) 6sin(2x) 2 4cos(2x) . /n ' 

; + - + i — '- + sin(2a;) 



/ 3 / 3 /■>•» 2 i. , 



f 

Jo 



(4.33) 



which has the window property 

" 30 dr 

— six) = 1 , (4.34) 
x 



the asymptotic behavior 

lim s(x) = -^x 5 + 0{x T ) , (4.35) 

and the inverse properties 

J x j *m *m) = - Q (4.36) 

or 

j o js{mwm = -e{i-k)--d{k-i) . (4.37) 

Then, the function s(x) allows us to invert Eq. (4.32) to get 

lnn 2 (£)= / s{k£) In P{k) . (4.38) 
Jo k 



4 See Appendix C for the equivalence between the single field reduction of Eq. (4.32) and the single field 
general slow-roll formula of Refs. [8-10] . 
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V. SUMMARY AND CONCLUSIONS 



We have developed a new 5N formalism for multi-component inflation that can be 
applied in the most general situations. We have introduced a perturbation variable SN 
whose form closely resembles the deviation of the e-folding number of expansion between 
two adjacent background trajectories in phase space, and whose evolution can be solved 
from sub-horizon to super- horizon scales. We have shown, though not rigorously yet, that 
SN(t- mi ) = <5jV(£fi n , t- mi ) on super-horizon scales, where 5N(t& n) t ini ) describes the perturba- 
tion of the e-folding number of expansion between a time i ini during inflation and the final 
comoving hypersurface t fin , on which the curvature perturbation 1Z C has achieved its final 
late time constant value. 

Taking the flat slicing during inflation, denoted by subscript f, we have derived a simple 
form for the evolution equation for SNf. This equation contains a source term that involves 
the scalar field perturbations explicitly. However, this source term is of first order in the 
deviations from scale invariance, and furthermore, its contribution to the spectrum of the 
curvature perturbation, which we have computed to leading order, can be expected to be of 
second order in the deviations from scale invariance. In the case when this source term can 
be neglected, which includes the case of conventional slow-roll inflation, we have obtained the 
power spectrum for 1Z C in a closed form in terms of functions determined by the background 
dynamics. In this case, we have found that the spectral formula is invertible, i.e. we have 
derived a formula that expresses the information of the inflationary dynamics directly in 
terms of the final curvature perturbation spectrum. 

The formula for the spectrum of 1Z C derived in this paper can be applied to a very wide 
class of models. Nevertheless, it will be useful to extend our formula to second order for 
greater accuracy and to include the possibility of super- horizon contributions [15]. Consider- 
ing special cases, such as the standard slow- roll limit, may help to make our formalism more 
accessible, and application to physically interesting models will also be useful. In addition, 
as mentioned above, we have not given a rigorous proof of the equivalence between 5N and 
5J\f on super-horizon scales. We plan to explore these issues in a forthcoming paper [13]. 
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APPENDIX A: CO VARIANT PARTIAL DERIVATIVES WITH RESPECT TO 

{<j>,4>) AND (0,tt) 

First we define the covariant partial derivatives in the phase space (0, 0) . For a scalar 
function F(0, 0) and phase space coordinates (0 a ,0 Q ), we have 

F = F^r + F^J a , (Al) 

where F^a and F^ a are the partial derivatives with respect to the coordinates <p a and a , 
and <j) a and <p a are the time derivatives of the coordinates 0" and Q . However, we wish to 
write all equations in a geometric manner. In order not to change the form of Eq. (Al), we 
should define the covariant partial derivatives as 

Fr = Fpe£ + F^<%, F^ = F^e a a , (A2) 

«V = V b e», (e^ b = 0, (A3) 

where the with a = 1, . . . , D are a complete set of basis vectors in the D-dimensional 
field space. Then we have 

F = F^ + F^ a , (A4) 
0; b = O, ^ = 55, (A5) 
Fcpjp — F^b = -F^ c i? c bda d , (A6) 

F^b^a. — F^b = F^ b< p a — -F^a^b = , (A7) 

where i? a b c 5 = (V a Vb — VbV a ) e 5 c is the curvature tensor on the scalar field space. 

Similarly, we may define the covariant partial derivatives in the phase space (0, n) , where 

^ = | • (A8) 

For a scalar function G(0, 7r), we define 

e a 

G^b. = G^e a + G n a— ^ , G^a = G^e^ , (A9) 



Then we have 



(e*V> = V b e», (e^ b =0. (A10) 

% = G^n a + G^ , (All) 

t#, = 0, Kb = 51, (A12) 

G^b^a — G^a^b = GTrci^bdaTT 11 , (^13) 

G-n-b^a — G^b = G^-b^a — G^a^b = . (A14) 
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The relationship between these partial derivatives is given by 



dQ_ 

d(j) a 



dQ_ 
<90 a 



and 



dQ_ 
<90 a 



In particular 
5Q = 



dQ_ 

dQ_ 

d<p a 

dQ_ 



+ 



+ 



+ 



where 



dQ 



dQ 



dQ 



dn 1 



1 (dQ , b 

7T 



H \d7i h 



4> 

dH 



dQ 



dQ 



d(j) a J w QH 2 \d7r h 



dQ \ (di: h 



d„ — 7T 



dir h 
1 f dQ 



H \d7T h 

1 ( dQ 



dH 



d<f) E 



H \d7i h 



rb 1 b 

a - -7T 7T a 



dQ 



3 H 2 \ dn a 



1 H 

1 + 



7T 



90' 

<97T a 



5vr a - 7T 8 
5vr a - 7T S 



if 3ff 



H 



5n a = 



H 



As noted in Section IIIB3, we may interpret 

5'<j) a - = 5 



dr 



(A15) 



(A16) 

(A17) 

(A18) 
(A19) 
(A20) 

(A21) 

(A22) 



which makes the form of SQ in Eq. (A18) identical to that of the perturbation of the back- 
ground quantity. Furthermore, noting that 



dH 

d^ 



+ 



dH 
<90 a 



-<f>*A)=HK + ^ j 1l, 



we find it is also possible to interpret 



5ir a - 7T 8 



A + K + ^K 



dAf 



(A23) 



(A24) 



at least on super-horizon scales or flat hyper surf aces. 
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APPENDIX B: PROPERTIES OF n (x), z\\(x), z ± (x), W(x) AND s(x) 



Here we list some mathematical properties of the functions used in Section IV. 
Definitions: 

1 



1 — ix — —x e 



n (x) 



, . 3sin(2a;) 6cos(2x) 4sin(2x) 

z \\ ( x ) = — dr^ zt^ z— 1 - + co < 2x ) ~ 1 



x° 



x z 



X 



, . 9 sin x 9 cos x 3 sin a; 

Z±(X) = g 

X 6 



X" 



X 



W(x) 



3sin(2:r) 3cos(2a;) 3sin(2x) 



2a; 3 



x/ 



2x 



- 1 



s(x) = - 



3 3cos(2:r) 6sin(2:r) 2 4cos(2x) , ' 

; + - + i — - + sin(2a;) 



/ y» 3 ,f> 3 / y* 2 , . i . 

tX-' JU JU JU JU 



Relations: 



n o n 'o + n = , 

x 



1 



Re [no(a;)] Im [n (a;)]' — Re [no(a;)]'lm [no(x)] = -x A 

9 



z\\(x) 



18 Im [n (x)} Re [n' (x)} = W{x) - - W'{x) 



x; 1 



9 



z±(x) = -^Im [n Q (x)} 



x 



Asymptotics: 



limno(x) = 1 + -x 2 + 0(x A ) + —x 5 + 0(ix 



x^O 



6 



45 



Kmz ll (x) = -x 2 + 0(x*) , 

Umz±(x) = \x 2 + 0(x 4 ) , 
x^o 5 

limW(a;) = -x 2 + 0{x A ) , 
4 

lim s(x) = x 5 + 0(x 7 ) . 

x^O 457T 



Window property: 



-z±(x) 



F 

Jo 

/ —[-xW\x)] = l, 
Jo x 

f°°dx 

/ — s(x) = l. 
Jo x 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 

(B6) 
(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 

(B15) 
(B16) 
(B17) 
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Degeneracy: 



f°° dx 1 , . r° dx 1 r , , . , 

/ Z \\{x)= [-xz|(a:)] = , 

Jo X X Jo X X 

r dx 1 un ^ 

/ W(x) = 

Jo XX 



[-xW'(x)\ = 



Inversion: 



/°° dk 
T s(ko Zll m = -9(t - o , 

/°° dk f 3 
T s(kQw(kZ) = -o(z-0-^o(c-z). 



(B18) 
(B19) 

(B20) 
(B21) 



where 6{x) = for x < and 0(a;) = 1 for x > 0. See Ref. [9] for a graph of the window 
function —xW'(x). 

APPENDIX C: CONNECTION WITH SINGLE FIELD FORMULA 

The single field general slow-roll formula for the curvature perturbation spectrum is [8-10] 

/°° dp r i i f 

-|[-^W'(^)][ln- + -^j , (CI) 

where 



/(lnO = 



H 



This is to be compared with the single field reduction of Eq. (4.32), 



InP 



o s 



, ^ 2n' 



where 



and 



n 2 



2 7 ra 2 C 2 / 
3 ^ rt dt 



a 



6H 



f - 



(C2) 



(C3) 



(C4) 



(C5) 



The single field formula can be inverted to obtain information about the inflaton dynamics 
from the primordial spectrum [11], 



This may be compared with the single field reduction of Eq. (4.38), 

f°° dk 

lnn 2 = / ^ s (k£,)\nP , 
Jo k 



(C6) 



(C7) 
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where m(x) and s(x) are related by 



m(x) — — m'(x) = s(x) + — s'(x) 



(C8) 



Thus the single field formulae and the single field reductions of the multi-field formulae 
are equivalent if 

2 17' 1 If 

+ + (C9) 

To check the above correspondence explicitly, we first note that in the single field case 
Eq. (3.13) reduces to 



HN: 



El 

A>2 



-)-■ 

6/ a 3 



Integrating this in time, we have 

Ni 



H 



C (El _ 1^1 



dt 



6/ a 3 



Therefore 



Ni = 



a 3 <j) 



H 



'ft'a 3 



1 

30 



1 + 



2cr0 



3^2 rt 



H 



H 2 dt 



H<j> ' j 02 



This gives 



lnNj + - 



2 W 



In I 

30 



3 ^ 

The other contributions to II are evaluated as 



2 _ 1 2 
3#0 ~ 3H 2 



In 



3// 
~2% 



2\ 2 



2(# 2 )' 

+ ~ 

3 H 2 



In 



3# 2 
~2V 



2\ 2 



2 0^ 
+ 3/J 2 



and 



Therefore 



In 



xaH 



2 

-2— . 

H 2 



, ^2 2n' , 

In IT + -— ~ In 



27T0/ 



2 5 2 
3^0 "3^ 



We see that this is identical to [9] 



, 1 2/' 



2 \ 2 



2tt0 



3 #0 3 H 2 



(CIO) 



(Cll) 



(C12) 
(C13) 

(C14) 

(C15) 

(C16) 
(C17) 



(C18) 
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APPENDIX D: CONFIGURATION SPACE FORMULA 



For /^-dimensional field space, the phase space is 2D-dimensional. On super-horizon scales 
during standard slow-roll inflation, however, half of the solutions are decaying, and eventually 
the phase space reduces to D-dimensions. Then is uniquely determined at each point in 
the configuration space, i.e. = 0(0). In such a case, we can define N((f>) = iV(0,0(0)) with 



V a iV = iV 0a + ^bV a b . 

Substituting into Eq. (3.30) gives 

6N = (V a iV) 50 a + iV; a \d'<p a - 0M - 50 b V b a 



Note that 



A a = 50 a - <j) a A - 



(Dl) 



:d2) 



(D3) 



is gauge invariant, and that it vanishes for the super- horizon adiabatic growing mode of 
Eq. (3.19). 

More generally, Eq. (3.2) gives 



3#V b a = -h^V^ - T^rV^ - V b a - ^0 a c V b c 



and we have 



V b a = -i? a cdb c d + 



(V b a )"+ (V b c ) (V c a ) . 
Therefore, using Eqs. (3.14) and (3.18), we have 



A a + 3#A a + ( V b a + ^0>b ) A b + q 2 S(Pt = . 



(D4) 
(D5) 

(D6) 



In the case of standard slow-roll, we see from Eqs. (D4) and (D5) that 

< 1 . 



(D7) 



Hence A a is decaying on super-horizon scales. Therefore the configuration space formula 
[3,6,7] 

5N = (V a iV) 50 a (D8) 
is valid for standard slow-roll on super-horizon scales. 
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